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SUMMARY 


A multiple-grid algorithm for use in efficiently obtaining steady solutions 
to the Euler and Navier-S tokes equations is presented. The convergence of a 
simple, explicit fine-grid solution procedure is accelerated on a sequence of 
successively coarser grids by a coarse-grid information propagation method which 
rapidly eliminates transients from the computational domain. This use of 
multiple-gridding to increase the convergence rate results in substantially 
reduced work requirements for the numerical solution of a wide range of flow 
problems. Computational results are presented for subsonic and transonic inviscid 
flows and for laminar and turbulent, attached and separated, subsonic viscous 
flows. Work reduction factors as large as eight, in comparison to the basic 
fine-grid algorithm, have been obtained. Possibilities for further performance 
improvement are discussed. 


INTRODUCTION 


When steady solutions to the Euler or Navier-Stokes equations are sought, it 
is common practice to compute them as temporal asymptotes to the unsteady 
equations of motion, subject to steady boundary conditions and the hypothesis that 
such solutions exist. This is done because the unsteady equations are either 
purely hyperbolic, in the case of the Euler equations, or hyperbolic - parabolic, 
for the Navier-Stokes equations, and are thus amenable to numerical solution by 
robust time-marching procedures. Such procedures are relatively easy to 
implement, readily capture weak solutions and, to the extent that the computation 
is time-accurate, allow a physical interpretation of the transient behavior, which 
may prove valuable in analyzing the performance of the algorithm. 

On the other hand, relaxation procedures for solving the steady versions of 
the Euler or Navier-Stokes equations are a topic of current research and have not 
as yet resulted in production algorithms. Attention has focused primarily on the 
steady Euler equations, where the research in progress includes the work of 
Bruneau, Chattot, Laminie and Guiu-Roux (1982) and Johnson (1982b). 

Both explicit and implicit time-marching procedures are presently in 
widespread use. Explicit methods, such as the MacCormack (1969) scheme, are 
simple, easily vectorizable and allow a good deal of flexibility in the treatment 




of boundary conditions. Their largest shortcoming lies in their conditional 
stability, which may place rather severe limitations on the time step size 
permissible on a given grid. When only a steady solution is sought and the 
accurate resolution of transients is of no consequence, explicit methods may 
consequently exhibit poor convergence rates. 

One remedy to the slow convergence of explicit schemes which has been 
investigated in recent years by Beam and Warming (1976), Steger (1977), MacCormack 
(1981), Lerat, Sid&s and Daru (1982) and others, is the construction of implicit 
methods. These methods, at the expense of a higher operations count, are 
generally unconditionally linearly stable and hopefully permit time steps to be 
taken as large as is consistent with accuracy requirements. In practice, large 
time steps may excite nonlinear instabilities and the choice of boundary condition 
implementation may introduce a stability limit. 

Another approach is to maintain the explicit nature of the integration scheme 
but to abandon the temporal accuracy requirement in favor of the introduction of 
convergence acceleration techniques, such as numerical methods with stability 
bounds in excess of the CFL limit or enhanced temporal damping procedures. Such 
techniques are exhibited in the work of Essers (1978), Viviand (1981), Jameson, 
Schmidt and Turkel (1981) and Rizzi and Eriksson (1981). 

Recently, efforts have been made at adapting multigrid techniques, which have 
shown substantial promise for the fast numerical solution of elliptic partial 
differential equations, to the Euler equations. Here, the nonlinear, nonelliptic 
nature of the equations, the existence of weak solutions and the fact that a 
coupled system of equations, rather than a single scalar equation, must be solved 
have been major obstacles. Jespersen (1983) has made notable progress in this 
direction. 

The present work, which is based on a concept introduced by Ni (1981), uses a 
sequence of successively coarser grids to accelerate the convergence of an 
explicit fine-grid solution procedure. In this general sense, it is a 
multiple-grid method. However, the details of its implementation are a good deal 
simpler than is the case with the conventional multigrid approach. The 
coarse-grid acceleration scheme presented here is quite modular. It may be used 
without modification in conjunction with a number of different fine-grid solution 
procedures and with any set of flow equations in the hierarchy ranging from the 
Euler equations to the full Navier-Stokes equations. 


EQUATIONS OF MOTION 


The nondimens ional equations of motion may be written in conservation-law 
form as 
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Here p, u, v, p, a and E are respectively density, velocity components in the x- 
and y- directions, pressure, sound speed and total energy per unit volume. The 
total energy per unit volume may be expressed as 

E = p(e + -|(u 2 + v 2 )) 

where the specific internal energy, e, is related to the pressure and density by 
the simple law of a calorically perfect gas 


P = (y - l)pe 


with y denoting the ratio of specific heats. The coefficient of thermal 
conductivity, k, and the viscosity coefficients, A and y, are assumed to be 
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functions only of temperature. Furthermore, by invoking Stokes’ assumption of 
zero bulk viscosity, X may be expressed in terms of the dynamic viscosity y as 



Re and Pr denote the Reynolds and Prandtl numbers, respectively. 

Although, for simplicity, the equations of motion are presented here written 
in Cartesian coordinates, Viviand (1974) has shown that their strong conservation 
law form may be maintained under an arbitrary time-dependent transformation of 
coordinates. Explicit detail concerning the generalized coordinate version of 
these equations which is employed in the computations to be discussed subsequently 
has been provided by Steger (1977) and need not be repeated here. 

We note that the thin-layer approximation, in the words of Baldwin and Lomax 
(1978), ”... evolves directly from a realistic assessment of what is really being 
computed in a typical high Reynolds number Navier-Stokes simulation.” A highly 
stretched mesh is used to resolve the large flow gradients normal to the 
vorticity-generating surface. Consequently, because of limitations on computer 
capacity, the diffusion terms involving derivatives parallel to the surface are 
not resolved well enough to merit their computation. 

Similar viscous terms are also neglected in the classical boundary layer 
approximation. However, while the boundary layer approximation replaces the 
normal momentum equation with the assumption that the normal pressure gradient is 
zero across the viscous layer, all momentum equations are retained in the 
thin-layer approximation and no assumptions are made concerning the pressure. 
Consequently, the separation point is not a singularity of the thin-layer model 
equations nor do the problems associated with matching a boundary layer solution 
to an inviscid outer flow occur when they are used. 

In practice, the thin-layer assumption is implemented by using a body-fitted 
coordinate system and neglecting the viscous terms in the coordinate direction 
along the body. For Cartesian coordinates, with x representing the 
body-conforming coordinate, the thin-layer version of the Navier-Stokes equations 
is as given above. 

The effects of turbulence are simulated by means of a two-layer algebraic 
eddy viscosity model. In the stress terms of the Navier-Stokes equations, the 
coefficient of dynamic viscosity, y, is replaced by y + y fc , where y t is the 
coefficient of eddy viscosity. Similarly, in the heat flux terms, the coefficient 
of thermal conductivity, k, is replaced by k + c p^t / Pr t> w ^ ere Pr t t ^ ie 
turbulent Prandtl number. The eddy viscosity is determined by the method of 
Baldwin and Lomax (1978) which is patterned after that of Cebeci (1970) with 
modifications to avoid the necessity of finding the boundary layer edge. Details 
may be found in the work of Baldwin and Lomax. 


FINE-GRID SOLUTION PROCEDURE 


The fine-grid integration scheme employed in this work is the two-step 
Lax-Wendroff method known as the MacCormack (1969) scheme. We observe that 
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schemes of Lax-Wendroff type may be arrived at intuitively by using Taylor's 
theorem to write the approximation: 

6q = Atq fc + q tt (2) 

where we define the "correction" to q such that 


6q = q (t + At) - q(t) 

Since we seek solutions to Eqn.(l), time derivatives may be expressed as space 
derivatives: 


q = -(F + G ) 
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where A and B are the Jacobian matrices: 
A = 3F/3q B = 3G/3q 


Substitution into Eqn. (2) results in: 


6q = -At(F x + G y ) 



One-Step Methods 

Second-order accurate spatial discretization of Eqn. (3) yields a one-step 
Lax-Wendroff method. Prior to discussing two-step schemes, we derive the one-step 
Lax-Wendroff fine-grid solution procedure used by Ni (1981). This is a necessary 
prerequisite to the development of the coarse-grid acceleration scheme. 

If we make the following finite-volume type approximations: 
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and define the "change" in q at cell centers such that: 
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it then follows that a discrete approximation to the first-order term in Eqn.(3) 
may be written as 
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Consistent with the above approximations and definition, we may write the 
approximation: 
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If we then approximate the second-order terms in Eqn. (3) such that: 
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we may combine Eqns. (5) and (6) to yield: 



Eqns.(4) and (7) constitute the one-step Lax-Wendroff method used as a basic 
integration scheme by Ni. He gives the following heuristic interpretation to 
these equations: the first calculates the change in q occurring in a control 
volume during the increment At, while the second distributes the effects of the 
changes occurring in four nearest-neighbor control volumes to their common central 
nodal point. There they are combined to form the correction to the vector of 
conservation variables, as illustrated in Fig. 1. This interpretation motivates 
the construction of the coarse-grid acceleration scheme to be discussed 
subsequently. 

Notice that Ni’s scheme may also be considered to be a two-step scheme with a 
full time increment predictor defined by Eqn.(4) and a corrector defined by 
Eqn.(7). However, such an interpretation would not be consistent with the general 
practice of avoiding the computation of Jacobian matrices in two-step schemes. 

Two-Step Methods 


Following Richtmyer (1962) many two-step Lax-Wendroff schemes have been 
developed. They have superseded the one-step schemes by virtue of their lower 
operations counts. MacCormack's method is a particularly popular and efficient 
member of this class. The forward predictor - backward corrector version of this 
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method may be written as 
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First derivatives in the viscous terms are backward differenced in the predictor 
and forward differenced in the corrector. 

The two-step methods of Burstein (1966) and Lapidus (1967) were also examined 
during the course of this work. Although they are spatially symmetric, their use 
resulted in no observable accuracy advantage over MacCormack* s method, for the 
test cases considered. Furthermore, they are computationally less efficient than 
MacCormack 1 s method. Details are contained in Johnson (1982a). 

This approach to solving fluid flow problems is quite robust and has been in 
widespread and successful use for some time, both for the time-accurate 
computation of unsteady flow and for the time-asymptotic solution of steady flow 
problems. In the latter case, where accurate resolution of physical transients is 
not required, the numerical stability limitation inherent in these explicit 
methods may severely restrict the speed of their convergence to the steady state. 
Providing a method to accelerate convergence in this case is the objective of this 
work. 


COARSE-GRID ACCELERATION SCHEME 


Given the fine-grid corrections, which may be computed by any one- or 
two-step Lax-Wendroff scheme, we wish to use successively coarser grids to 
propagate these corrections throughout the computational domain, thus accelerating 
convergence to the steady state while maintaining the accuracy determined by the 
fine-grid discretization. Given a basic fine grid with the number of points in 
each direction expressible as n(2 p ) + 1 for p and n integers such that p « 0 and 
n > 2, where p is the number of grid coarsenings and n is the number of 
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coarsest-grid intervals, let successively coarser grids be defined by successive 
deletion of every other point in each coordinate direction. 

Full Coarse-Grid Scheme 


The full coarse-grid acceleration scheme, as illustrated in Fig. 2, replaces 
the computation of coarse grid changes by Eqn. (4) with a restriction of the latest 
fine-grid correction. This restricted fine-grid correction is then distributed 
according to a coarse-grid version of Eqn. (7) to obtain a coarse-grid correction. 
This coarse-grid correction is, in turn, prolonged to the fine grid to become the 
new fine-grid correction. One time-cycle of the multiple-grid scheme is composed 
of an application of some Lax-Wendroff method on the fine grid followed by an 
application of the coarse-grid propagation scheme to each successively coarser 
grid. The flow of information in this process is depicted in Fig. 3. 

In the basic integration scheme, a change at one grid point affects only its 
nearest neighbors while, in a k-level multiple-grid scheme, the same change 
affects all points up to 2 k-1 mesh spacings distant. Furthermore, since the 
change is always determined by information from the fine grid and simply 
propagated by the distribution formulae for coarser grids, fine grid accuracy is 
maintained. 

This concept for convergence acceleration was introduced by Ni for use in 
conjunction with his one-step Lax-Wendroff scheme, as described above. He 
illustrated its utility by solving the homoenthalpic two-dimensional Euler 
equations. Johnson (1982a) extended its applicability to the class of two-step 
Lax-Wendroff schemes and solved the full two-dimensional Euler equations. 

Convective Coarse-Grid Scheme 


In Johnson (1982c), consideration of the physical processes being modelled in 
a viscous flow computation led to the formulation of an alternative coarse-grid 
scheme. Dissipative effects have a local character and their influence need not 
be taken into account in the construction of coarse-grid distribution formulae. 
Rather, it is the convective terms, with their global character, which are the key 
element in coarse-grid propagation. Hence, a coarse-grid scheme for viscous flow 
computations may be formulated on the basis of the inviscid equations of motion. 
Such a convective coarse-grid scheme is inherently more efficient than the full 
coarse-grid scheme because of the diminished computational effort associated with 
forming the Jacobian matrices of the Euler flux vectors rather than those of the 
viscous flux vectors. An additional benefit is that the convective coarse-grid 
scheme leads to a multiple-grid convergence acceleration procedure which is 
independent of the nature of the dissipative terms retained in the viscous model 
equations. That is to say: the coarse-grid scheme based on the Euler equations 
may be employed, without modification, to accelerate the convergence of viscous 
flow computations based on the Navier-Stokes equations, the thin-layer equations, 
or any other viscous model equations which contain the full inviscid Euler 
equations . 

The correctness of the heuristic physical reasoning used in the formulation 
of the convective coarse-grid scheme is supported by the computational results for 
viscous flow to be presented subsequently. 
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COMPUTATIONAL RESULTS 


To investigate the robustness of the coarse-grid acceleration scheme and to 
explore the generality of its applicability, a large number of computational 
experiments have been performed. A representative sampling of those computations 
is reported here. 

The full Euler equations have been solved for both subsonic and transonic 
flow. The thin-layer version of the Navier-Stokes equations has been solved for 
attached and separated, laminar and turbulent, subsonic flows. All computations 
have been performed in two dimensions. Extension of the technique to the full 
Navier-Stokes equations or to three dimensions presents no essential difficulties 
and should, in fact, yield greater advantage than is illustrated by the present 
results. 

Problem Specification 

We consider the inviscid flow through an unstaggered cascade of bicircular 
arc airfoils at zero angle of attack, as illustrated in Fig. 4, and the viscous 
flow through a similar cascade of sting-mounted airfoils, as shown in Fig. 5. The 
boundary conditions used are also indicated in these Figures. At the upstream 
domain boundary, total pressure, total temperature and flow angle are specified. 

At the downstream boundary, the static pressure is fixed. Along inviscid lateral 
boundaries, the tangency condition is applied, while, along solid walls, the 
no-slip condition is applied and the temperature specified. The values of the 
dependent variables on the domain boundaries are updated only during the fine-grid 
computations. This decouples the coarse-grid acceleration scheme from the details 
of boundary condition implementation. Uniform flow at the isentropic Mach number 
implied by the ratio of exit static pressure to upstream total pressure is used as 
an initial. ;state. 

For the subcritical cases, the ratio of exit static pressure to upstream 

total pressure is 0.8430191, yielding an isentropic upstream Mach number of 0.500, 
while, for the supercritical cases a ratio of 0.7369520 is employed, resulting in 
an upstream Mach number of 0.675. In the viscous cases, the Reynolds numbers, 
based on cascade gap and critical speed, span the approximate range from 8.4 x 10 3 
to 2.0 x 10 5 . 

Symmetry and periodicity are invoked to limit the size of the computational 
domain. All inviscid computations were carried out using the 65 x 17 node fine 
grid as shown in Fig. 6 and the coarser grids as indicated in Table I. As is also 
illustrated in Fig. 6, three different fine grids were employed in the viscous 
flow computations. All have 65 x 33 nodal points and have their transverse grid 
lines located at the same positions. They differ in the positioning of their 
lateral grid lines. These are smoothly stretched away from the solid boundary in 
a geometric progression, starting from the initial spacings indicated in Fig. 6. 
These fine grids each allow the construction of four successively coarser grids, 
as shown in Table II. 
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Inviscid Flow 


The results of typical subcritical inviscid flow computations are illustrated 
in Fig. 7. The airfoils have a thickness-to-chord ratio of 0.200. Isomach 
contours and the surface Mach number distributions along the upper and lower 
domain boundaries are shown. We observe that the fore-and-aft geometric symmetry 
is well preserved in the computed flow field. 

All computations reported here are considered to be converged when the _ 6 
average absolute value of the unsealed fine-grid residual in pu drops below 10 , 

representing a decrease of approximately five decades from its original value. 

This convergence criterion is more than sufficient for engineering purposes and 
could be relaxed by at least two decades without any significant accuracy 
deterioration in the solution. The converged single- and multiple-grid results 
were, as should be expected, identical, confirming that the solution obtained is 
not a function of grid sequence length. 

Based on the computational work required to produce a converged solution, we 
may estimate a multiple-grid work reduction factor and a corresponding optimal 
grid sequence length. We define the work reduction factor to be the ratio of the 
work required to produce a converged solution using a single fine grid, to the 
work required to produce the same result using that multiple-grid sequence length 
which minimizes the computational work. The results obtained are recorded in 
Table III. For the case at hand, a sequence of three grids produced a work 
reduction factor of 5.1. The single- and optimal multiple-grid convergence 
histories are also shown in Fig. 7. 

Supercritical results for the same thickness-to-chord ratio are shown in 
Fig. 8. Here, we observe that the shock is crisply resolved by the multiple-grid 
scheme although the explicit artificial viscosity term added to the algorithm for 
this purpose is only applied on the fine grid. We hasten to add that no artificial 
viscosity beyond that inherent in the basic multiple-grid algorithm was employed 
in either the inviscid or viscous subcritical computations reported here. A 
sequence of four grids produced a work reduction factor of 2.1 for this inviscid 
supercritical flow case. 

Viscous Flow 


Viscous flow computations have been performed for the combinations of 
Reynolds number, thickness-to-chord ratio and fine grid configuration indicated in 
Table III. A sampling of the results is illustrated in Figs. 9 through 12. In 
each case, isomach contours, u-velocity profiles and the single- and optimal 
multiple-grid convergence histories are shown. Note that the contour levels 
displayed are not equally spaced and that the same levels were used in all 
illustrations. Nevertheless, they provide a good qualitative indication of the 
nature of the computed flowfields. More quantitative information is provided by 
the normalized u-velocity profiles. The u-velocity, normalized with its value at 
the domain top boundary and same streamwise station, is plotted as a function of 
relative distance from the airfoil chord line. Curves for every second streamwise 
station, starting with the airfoil leading edge and ending at the outflow 
boundary, are displayed. They are staggered in proportion to the spacing of their 
respective streamwise stations. 
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The thickness-to-chord ratios employed ranged from 0.100 to 0.000. Both 
laminar and turbulent flows have been computed. In the turbulent cases, the flow 
was tripped at the airfoil leading edge. Attached and separated cases were 
computed for both laminar and turbulent flow. 

The work reduction factors obtained ranged from 1.7 to 8.2 and are recorded 
in Table III, together with the optimal grid sequence lengths. We note that 
Johnson (1982c) observed that although the work reduction factor and, possibly, 
the optimal grid sequence length decrease with increasing grid stretching, they do 
not appear to decrease with increasing Reynolds number. 


ALGORITHMIC CONSIDERATIONS 


During the course of this work, several issues bearing on the speed, 
computational efficiency, storage requirements and general applicability of this 
multiple-grid convergence acceleration scheme have come to mind. As most of these 
issues deserve a good deal more study, we let the following brief remarks suffice 
for the present. 

Boundary Condition Implementation 

As mentioned previously, boundary conditions are only enforced on the fine 
grid. This has the advantage of decoupling the coarse-grid acceleration scheme 
from both the physical and numerical nature of these boundary conditions. That is 
to say: the coarse-grid scheme always sees a Dirichlet problem. While this 
enhances the modularity of the coarse-grid scheme, it is not clear that such a 
boundary condition treatment serves to maximize the convergence acceleration of 
the resultant algorithm. 

Accelerator Linearization 


The computation of the Jacobian matrices used in the coarse-grid acceleration 
scheme has a non-trivial influence on the coarse-grid operations count and, hence, 
on the work reduction factor. We have found that the Jacobians need not be 
updated at each time cycle, but may be lagged by a substantial margin without any 
adverse effect on the resultant convergence history. In fact, in some cases this 
improves the convergence rate and, in so doing, raises the possibility of 
tailoring the Jacobian matrices to accelerate convergence rather than calculating 
them based on the actual state of the flow being computed. 

Storage Minimization 

Linearizing the coarse-grid acceleration scheme in the sense described above 
implies storing the Jacobian matrices rather than recomputing them based on the 
latest available flow state. Should this storage requirement constitute an 
unacceptable burden, it may be essentially eliminated by periodically storing 
information on the flow state at a small number of representative points and 
reconstructing a set of Jacobians from this information at each time cycle. Such 
a strategy would be in complete accord with our linearization ideas. 
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Restriction and Prolongation Operators 


In the present computations, injection is used as the restriction operator 
and linear interpolation is chosen as the prolongation operator. These choices 
may not be optimal for use on highly stretched grids. Better choices could 
increase both the optimal grid sequence length and the work reduction factor. 

Coarse-Grid Formation Strategy 

Similar consequences might result from an alternative coarse-grid formation 
strategy. Given the convective nature of the coarse-grid scheme used here, it is 
entirely plausible that, particularly in viscous flow computations, the coarser 
grids should have a more uniform structure than is obtained by successive deletion 
of every other grid line from a highly stretched fine grid. 

Coarse-Grid Propagation Strategy 

In distinction to conventional multi-gridding, the current algorithm cycles 
from the finest grid to the coarsest at each step to the solution procedure. 
Consequently, relatively many computations are performed on the finest grid. 
Computational expense could be reduced by developing an alternative coarse-grid 
propagation strategy which would spend more of its time at the coarser grid 
levels. 

Vectorization 


As both the fine-grid solution procedure and the coarse-grid acceleration 
scheme used here are explicit, the resultant multiple-grid algorithm may be 
readily vectorized. Such vectorization is presently being performed and will be 
reported at a later date by Chima and Johnson (1983). 

Applicability of the Accelerator 

As the coarse-grid scheme uses a one-step Lax-Wendroff framework and as this 
concept was first employed by Ni to accelerate his one-step Lax-Wendroff fine-grid 
solution procedure, it is natural that the first generalization should be to the 
class of two-step Lax-Wendroff schemes. However, there is no compelling reason to 
expect that its applicability is limited to this class of methods. We expect that 
it may be used with a variety of explicit fine-grid solution procedures and, in 
fact, Stubbs (1983) has demonstrated its use with the implicit MacCormack (1981) 
scheme. 


CONCLUSIONS 


An explicit multiple-grid algorithm for the efficient solution of the steady 
Euler and Navier-Stokes equations has been presented. 

The coarse-grid scheme used to accelerate convergence is compatible with all 
one- and two-step Lax-Wendroff algorithms. Here, it has been used in conjunction 
with MacCormack’ s method. 
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The convective version of the coarse-grid scheme may be used with any set of 
flow equations in the hierarchy ranging from the Euler equations to the full 
Navier-Stokes equations. Here, we have used it with the thin- layer version of the 
Navier-Stokes equations. 

Computational results have been presented for subsonic and transonic inviscid 
flows and for separated and attached, laminar and turbulent, subsonic viscous 
flows. 

Multiple-gridding yielded work reduction factors ranging approximately from 
two to eight, over a fairly broad range of flow conditions. Ample margin remains 
for further performance improvements. 

Given the encouraging results obtained to date, more comprehensive testing 
and more sophisticated applications of the inviscid and viscous flow convergence 
acceleration ideas presented here are planned. 
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TABLE I 


INVISCID FLOW GRID DESCRIPTIONS 


Grid 

i 

2 

3 

A 

Number of 
Points 

65 x 17 

33 x 9 

17 x 5 

9x3 


TABLE II. - VISCOUS FLOW GRID DESCRIPTIONS 


Grid 
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Number of 
Points 

65 x 33 
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17 x 9 

9x5 

5x3 






















TABLE III 


SUMMARY OF RESULTS 



Comments 


Inviscid 
Subcritical 
Fig. 7 


Re - 8.4*10 3 
Laminar 
Fig. 9 


Re « 8.4*10 3 
Turbulent 


Re - 3.4x10** 
Laminar 


Re - 2.0x10 s 
Turbulent 


Re - 8.4xl0 3 
Laminar 


Re - 3.4x10** 
Laminar 
Separated 
Fig. 10 


Re - 3.4x10** 
Turbulent 
Separated 
Fig. 11 


Re - 3.4X10** 
Laminar 


Re - 2.0x10 s 
Turbulent 
Fig. 12 


Thickness 

Initial 

Transverse 

Spacing 

Chord 

0.200 

0.03125 

0.200 

0.03125 

0.000 

0.01250 

0.000 

0.01250 

0.000 

0.00625 

0.000 

0.00250 

0.100 

0.01250 

0.100 

0.00250 

0.100 

0.00250 

0.050 

0.00625 

0.025 

0.00250 


Optimal 

Sequence 

Length 


Work 

Reduction 

Factor 
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FIGURE 1. - One-Step 

Lax-Wendrof f Scheme 


FIGURE 2, - Coarse Grid Scheme 
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FIGURE 3. - Multiple-Grid Information Flow 
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Inviscid Bicircular Arc 
Cascade Problem 





FIGURE 5. - Viscous Bicircular Arc 
Cascade Problem 
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Fine Grid Configurations 





Norm of Residual 


Isomacns 



Normalized u-velocity 
Velocity Profiles 


Convergence Histories : 



5000 

10000 

0 

5000 

Time Cycles 

D 

scaled p 

Time Cycles 


o 

scaled pu 



A 

scaled pv 



+ 

scaled E 



X 

average unsealed pu 
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FIGURE 10. - Viscous Flow, Re = 3.4 x 10 4 , Laminar, Separated 
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FIGURE 11. - Viscous Flow, Re = 3.4 x 10 4 , 
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